Sub-milliarcsecond Imaging of Sgr A* and M 87 
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Abstract. We present and discuss new result from mm-VLBI observations of M 87 and Sgr A*. 
The imaging of these sources with a spatial resolutions of a few to a few ten Schwarzschild radii 
offers new possibilities to study the immediate environment of super-massive black holes. 

1. Introduction 

In astronomy, the highest angular resolution is obtained with Very Long Baseline Interferometry 
(VLBI). VLBI observations at short millimeter wavelengths (mm-VLBI) provide angular 
resolutions in the ten micro-arcsecond range (i.e. 14 /nas for a 10, 000 km baseline at A = 1.3 mm 
[y = 230 GHz)). While VLBI observations at 2 mm and 1 mm wavelength are still in a stage of 
technical development, 3 mm VLBI experiments are performed on a regular basis. 

The Global mm-VLBI Array (GMVA 1 ) is operational since early 2000. At 86 GHz (A = 
3.5mm), it combines the European antennas (Effelsberg 100m, Pico Veleta 30m, the phased 
Plateau de Bure interferometer 6x15 m, Onsala 20 m, Metsahovi 14 m) with the VLBA. With 
the participation of the two sensitive IRAM telescopes and the 100 m Effelsberg telescope, the 
array sensitivity is improved by a factor of 3 — 4, when compared to the VLBA alone. For 
compact galactic and extragalactic radio sources, the GMVA provides VLBI images with an 
angular resolution of down to 40 //as. 

In the following we present and discuss results from mm-VLBI observations of the nuclei of 
M87 and of Sgr A*. When expressed in terms of Schwarzschild radii (R$ = 2GMbh/c 2 ), the 
spatial resolution obtained for M87 and Sgr A* become very comparable. M87 is about ~ 10 3 
times more distant than Sgr A*, but has a putative black hole at its center, which is also ~ 10 3 
times more massive than that of Sgr A*. Present day 3 mm-VLBI therefore allows to image 
both sources with a spatial resolution of only 10 — 20Rs- In view of the morphological difference 
between Sgr A* and M87 (the latter exhibits a long relativistic jet), a comparison of the sub- 
milliarcsecond structure and kinematics of both objects may lead to a better understanding of 
the astro-physical processes acting in the vicinity of such super-massive black holes. 



2. 3 mm VLBI observations of M 87 

The giant elliptical galaxy M 87 (NGC 4486, 3C 274, Virgo A) is one of the closest radio galaxies 
with a prominent radio jet. Its relative brightness at millimeter wavelengths makes it ideally 



1 web link: http://www.mpifr-bonn.mpg.de/old_mpifr/div/vlbi/globalmm 



suited for high angular resolution mm-VLBI studies, which address the still unanswered question 
how the powerful jets in radio-galaxies and quasars are launched and collimated. 

At a distance of D=16.75 Mpc (Whitmore et al. 199) an angular scale of 0.1 mas corresponds 
to a spatial scale of 2.5- 10 16 cm (0.0081 pc). This translates to 27.8Rs Schwarzschild radii for an 
assumed 3 • 10 9 M Q super- massive black hole at the center of M87 (e.g. Machetto et al. 1997). 
Global 3 mm VLBI observations, performed at the edge of the accessible VLBI observing bands, 
have revealed new images of the inner-most region of M87. The 3 mm VLBI images shown 
in Figure 1 were made from observations performed during 2002-2004 and show a one-sided 
core-jet structure, extending ~ 2 mas to the west. The complexity of the source structure and 
the only moderate dynamic range of the maps (lower sensitivity in earlier observations) do not 
yet allow to relate the individual structural components seen at the different epochs with each 
other and by this determine the jet kinematics. Despite the obvious limitations, which will be 
overcome in the near future, a number of basic statements can already be made: 

(i) In the maps, an unresolved compact component at the eastern end of the jet is seen. An 
upper limit to its size is given by the (elliptical) VLBI observing beam. The uniformly weighted 
and un-tapered beam size of the observations of April 2004 is (193 x 54)^as, which gives an 
upper limit to the size of the jet base (VLBI core) of (54 x 15)Rs- Circular Gaussian model 
fitting yields a typical flux density of ~ 0.5 Jy and a size of < 50fias (< 14i?s) for the VLBI 
core. This results in a brightness temperature of > 10 10 ' 5 K, clearly indicative of non-thermal 
(synchrotron) radiation. 

(ii) The 86 GHz VLBI map of April 2004 looks remarkably similar to two previously observed 
43 GHz VLBI maps, which indicate conical jet-expansion on sub-milliarcsecond scales (Junor et 
al. 1999, Ly et al. 2004). The new 86 GHz maps, which have a ~ 2 times higher resolution, 
reveal a bifurcated (possibly edge-brightened) expanding one-sided jet, with the southern edge 
of the emission cone being brighter and more compact than the northern part. The fact that 
the northern edge of the cone (at r < 0.5 mas) is almost invisible in the maps of 2002 and 2003 
may be explained with dynamic range limitations in these two maps. 

(iii) The large morphological differences seen between the three 86 GHz maps indicate structural 
variations on time-scales of less than one year. On angular scales < 0.5 mas, all three maps show 
some south-west oriented jet emission oriented along a position angle of ~ —110°. In April 2003, 
the source was ~ 1.5 — 2 times brighter than usual, and a bright component was seen at the 
map center (at zero position). Since in hybrid VLBI imaging using closure phases the absolute 
position information is lost, this bright component most likely should not be identified with the 
jet base (VLBI-core). Instead, the fainter component located ~ 0.2 mas in the north-east, may 
be identified with the VLBI-core. We note the tentative detection of a possible counter-jet at 
43 GHz (Ly et al. 2004), which however was not seen earlier (Junor et al. 1999), nor at longer 
wavelengths (Dodson et al. 2006), nor in our 86 GHz VLBI maps of 2002 and 2004. 

(iv) Owing to the limited time sampling at 86 GHz, a detailed kinematic study of the component 
motion in the inner jet of M87 is not yet possible. Similar to the case of many other VLBI 
jets, we may assume that the bright central component seen in April 2003 was ejected from 
the VLBI core after a flux density outburst. This would give an apparent rate of change of 
^ 0.2mas/yr, corresponding to an apparent velocity of > 0.05 c. This speed would be consistent 
with the observed motion of 0.04 ± 0.02 c seen at larger core separations with VLBI at 15 GHz 
(Kellermann et al. 2004), and with a lower limit of 0.25 c from VLBI at 43 GHz (Ly et al. 2004b). 
On the other hand, if we compare the relative positions of the two main secondary components 
seen in April 2002 and April 2004 (at r ~ 0.25 mas &: 0.50 mas, respectively) stationarity of the 
inner jet components cannot be excluded. More frequent VLBI observations are necessary to 
determine the jet kinematics. 

We like to point out that a measurement of VLBI core size and of the jet's transverse width 
as a function of core separation, are very important parameters for theoretical jet models. The 



class of models, which are based on the magnetic sling shot mechanism (Blandford & Payne 
1982) predict a minimum transverse width of the jet after initial acceleration to be of order 
> j-2Rl, (with bulk Lorentz factor 7 and light cylinder radius Rl, e.g. Tomimatsu & Takahashi 
2003). Notably such jets are expected to be hollow (edge-brightened), just as observed here. 
After collimation, the light cylinder has a radius of at least > (10 — 50) Rs (Camenzind 1990, 
Fendt & Memola, 2001, Casse & Keppens 2004), which sets a lower bound to the observable 
jet width. With jet Lorentz-factors in the range 3 < 7 < 6 (Biretta et al. 1999), we expect a 
minimum jet width of > (60 — 600)Rs- The observed diameter of the jet (at 0.5 < r < 2 mas) is 
2R L = (0.5 - 0.7) mas, or (140 - 200)-Rs (see Figure 1). This is at the lower end of the expected 
range, favoring a moderately low bulk Lorentz- fact or (7 < 2 — 3) in this region. 

At its origin near the black hole, the diameter of the jet depends on the size of the region, 
where it is initially accelerated. In the Blandford & Payne type models, the twisted magnetic 
field lines are anchored in the inner part of the rotating accretion disk. The minimum jet width 
is then determined by the diameter of the least stable orbit, which is at 6Rs for a Schwarzschild 
BH, and 2R$ for a rotating Kerr BH (e.g. Meier et al. 2001). The observed small size of the 
VLBI core (< 15-Rs) is already within a factor of 2 — 3 near these values, suggesting that the 
width of the 'jet-nozzle' is indeed of order of the light cylinder diameter. Since the radiating 
jet boundary cannot be indefinitely thin, one could expect that the observed jet width is a bit 
larger than the light cylinder. In this case, the observed small size of the VLBI core would point 
more towards jet models, in which the jet gains energy directly from BH-rotation (Blandford 
&: Znajek 1977, see also McKinney 2006). Better VLBI images and more specific theoretical 
modeling of jet opening angle and speed will help to discriminate between the different types 
of models (e.g. Blandford & Payne versus Blandford k, Znajek), and may tell to what degree 
black hole rotation plays a role in jet formation. 




Figure 1. 86 GHz VLBI maps of the jet in M87. The observations were performed (from left 
to right) in April 2002, 2003, and 2004. The maps result from global VLBI observations using 
the GMVA, which combines European radio telescopes with the VLBA. 



3. VLBI observations of Sgr A* at 3, 2, and 1 mm wavelength 

3 mm: From the early 1990's onwards, Sgr A* was repeatedly observed with VLBI at 3.5 mm 
(86 GHz) in various array configurations, using European, American and global VLBI arrays 
(Krichbaum et al. 1994, Rogers et al. 1994, Doeleman et al. 2001, etc.). In Figure 2 we show one 
of the early maps from global 3 mm VLBI, which reveals a partially resolved point sources, similar 
to the recently published VLBA image (Shen et al. 2005). At mm- wavelengths, the accurate 
determination of flux density and size is often affected by calibration uncertainties (low source 
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Figure 2. 86 GHz VLBI image of SgrA* 
observed in April 1997. The participating 
antennas are: Effelsberg (100 m), Pico Veleta 
(30 m), Haystack (37 m), Quabbin (14 m), Pie 
Town (25 m), Kitt Peak (12 m), Owens Valley 
(5 x 10 m). A circular Gaussian parameterizes 
SgrA* with flux density S = 1.73 ± 0.25 Jy 
and size 9 = 0.28 ± 0.08 mas. 
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Figure 3. Closure phase of SgrA* plotted 
versus time. Data are from a 3 station 
VLBI experiment performed in October 2005 
at 86 GHz using the telescopes at Effelsberg, 
Pico Veleta and Plateau de Bure (phased). 
The closure phase is zero within ±10°, 
suggesting a point-like or almost perfectly 
symmetrical brightness distribution. 
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Figure 4. Angular size of 
SgrA* plotted versus wave- 
length. Lines denote fits 
(aA^), excluding the short 
wavelengths A < 7 mm, which 
are affected by interstellar 
scattering. The fit param- 
eters a, (3 are given inside 
the figure. Published 3 mm 
measurements are plotted as 
open squares, new 3 mm data 
(this paper) by open circles. 
New measurements at 2 mm 
(shaded triangle, an upper 
limit) and 1.4 mm (diamond) 
are also shown. 



elevation on the northern hemisphere, atmospheric opacity variations). This requires special 
observing and calibration methods. In parallel to the ongoing research with the VLBA, we have 
continued to observed SgrA*, using the sensitive and large antennas at MPIfR (Effelsberg 100 m 
telescope) and at IRAM (Pico Veleta 30m, Plateau de Bure 6 x 15m). With baseline lengths 
of ~ 750 — 1700 km, and a detection sensitivity of > 50mJy (7a, 512 Mbit/sec) per baseline, 



this array measures the visibility amplitudes and (closure-) phases with particular high SNR. 
In April and October 2004, Sgr A* was moderately bright with flux densities of 2.40 ± 0.45 Jy 
and 2.10 ± 0.23 Jy, respectively. In October 2005, however, the source was detected with a much 
higher signal-to- noise ratio than before (SNR = 96 on PV-PdB). A preliminary analysis of the 
new data yields (via circular Gaussian model fitting) a flux density of S = 3.6 ± 0.5 Jy and 
a size of 9 = 0.27 ± 0.05 mas. The relatively high source flux is confirmed independently by 
ATCA measurements, performed 3 days after our VLBI experiment (S = 2.0 — 3.5 Jy, Z.-Q. 
Shen, priv. comm.). The source size appears marginally larger (by 2 — 3<j), when compared 
to our previous size measurements in 2004 (9 = 0.19 ± 0.05 mas and 0.18 ± 0.01 mas), and in 
comparison to the VLBA measurements of Shen et al. 2005, who obtained 9 = 0.18 ± 0.02 mas 
with the VLBA. As in previous observations, the closure phase (see Figure 3) is close to zero 
(within ±10°), showing no evidence for a deviation from the known point-like (or symmetric) 
structure. In Figure 5, we plot all available 86 GHz VLBI size measurements, using the data 
from the literature and including own unpublished old and new data. If one excludes the new 
data point of October 2005, the weighted mean average size is 0.18 ± 0.01 mas. This appears 
significantly smaller than the October 2005 measurement. For comparison we also show in Fig. 
5 the sizes at 43 GHz (weighted mean: 0.55 ±0.03 mas) and the flux density of Sgr A* at 86 GHz. 
It will need a denser 3 mm VLBI monitoring to check, if the brightness distribution of Sgr A* 
indeed is variable. In view of the observed flux density variations at millimeter wavelengths (see 
Miyoshi et al., Miyazaki et al., Marrone et al., this conference), structural variability of Sgr A* 
would be not unexpected (see Bower et al. 2004 for possible size variations at 7mm). 

2 mm: VLBI observations at wavelengths shorter than 3 mm are still difficult and are limited 
by the number of available radio telescopes and their sensitivity. In a pilot VLBI experiment 
at 147 GHz (2 mm), which aimed at first fringe detection across the Atlantic at this short 
wavelength, Sgr A* was observed among a dozen other AGN (Krichbaum et al. 2002). Several 
AGN were detected on the long (~ 8500 km) baseline between Spain (Pico Veleta) and the 
Arizona telescopes (Kitt Peak (KP) 12 m, Heinrich-Hertz telescope (HHT) 10 m), which limits 
the size of their compact emission regions to < 25 — 40 fias. Sgr A* was detected with a signal- 
to-noise ratio of ~ 7 on the short (85 MA) baseline between KP and HHT. With a total flux 
density of ~ 1.9 Jy and a correlated flux density of S c = 1.6 ± 0.25 Jy an upper limit to its size 
of < 0.7 mas can be derived, fully consistent with the measured sizes at 3 mm and 1.4 mm. 

1 mm: In April 2003, our attempt to observe Sgr A* with global VLBI at 230 GHz failed due 
to weather and technical problems. In this experiment, which was motivated in a similar manner 
as the 2 mm VLBI experiment described before, several bright AGN were observed. Finally two 
blazars (3C 454.3, 0716+714) were detected with SNR Z 7 on the 6.4 GA baseline between Pico 
Veleta and HHT (Krichbaum et al. 2004). Although limited by receiver performance and 
stability of local oscillators at HHT, this experiment demonstrates the technical feasibility of 
global VLBI at 1.3 mm, and sets a new record in angular resolution. For the detected compact 
radio sources we obtained sizes of < 30 fias. With such a high angular resolution, which would 
correspond to a spatial resolution of ~ 3 Schwarzschild radii in Sgr A*, the direct (VLBI-) 
imaging of the immediate vicinity of a SMBH should become possible within the next few years. 

We now like to make a remark with regard to the previous 215 GHz VLBI observations of 
Sgr A* published earlier (Krichbaum et al. 1998). Recently, Eckart et al. 2006 (also Eckart et 
al., this conference) detect an emission component Dl of ~ 170 mas size, located 80 ± 50 mas 
west of Sgr A*. With single dish radio telescopes operating at short mm-/sub-mm wavelengths, 
in-beam confusion between Sgr A* and Dl leads to an over-estimate of the total flux density 
of Sgr A* in this spectral band (by AS215GHZ — 0.75 Jy for T ~ 10 3 K), at least partially 
removing the so called sub-mm excess (e.g.. Serabyn et al. 1997, Falcke et al. 1998). In fact, a 
less pronounced or non-existing sub-mm excess, leads to a flat synchrotron-self-absorbed radio 
spectrum as observed in many other AGN. This then would allow an interpretation of the high 
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Figure 5. Measured VLBI sizes of Sgr A* at 
86 GHz (filled circles) and at 43 GHz (open 
diamonds) plotted versus time. Data are 
taken from the literature, including our new 
and unpublished 86 GHz data. All sizes were 
calculated from the geometric mean of major 
and minor axis. We also show the VLBI flux 
density (filled squares) at 86 GHz. 



Figure 6. Right: Simulated visibility amplitudes (top) and closure phase (bottom) for a future 
1.3 mm VLBI experiment on Sgr A* with the following telescopes: Pico Veleta (Spain), Plateau 
de Bure (France), CARMA (California), SMA (Hawaii), HHT (Arizona) and APEX (Chile). 
For the simulation, a circular Gaussian of 2.5 Jy and ~ 25 /xas FWHM size was assumed. 




energy data via the synchrotron-self-Compton process (see Eckart et al. 2004 & 2006). 

The existence of such an extended emission region, which is fully resolved and therefore 
invisible to VLBI, was considered by Krichbaum et al. 1998 as a possibility to explain the 
apparent discrepancy between observed total flux (from single antenna measurements) and 
the correlated VLBI flux, which originates from the compact component. Relating the sub- 
mm excess of the total spectrum mainly to the flux contribution of component Dl, allows 
to redetermine the total flux of Sgr A* at 215 GHz and reject 'case 1' in the discussion of 
Krichbaum et al. 1998. This leads to refined estimate of the observed size for Sgr A* at 215 GHz: 
9 = (110 ± 60)^as, which corresponds to a linear size of Rs = 11 ± 6 Schwarzschild radii and a 
robust lower limit of Rs > 5. 

4. Intrinsic Source Size 

The observed frequency dependence of the source size of Sgr A* (Fig. 4) is commonly interpreted 
as being due to scatter broadening by the intervening interstellar medium. The intrinsic size 
is determined by the deconvolution of the observed size # bs with the extrapolated scattering 

size 6 scat - #int = \J #obs ~~ ^scat • The exact form of the assumed scattering law (0 sca t = aX 13 , 
with fit parameters a, f3) determines #i n t. In Figure 4 we show a few illustrative least-square fit 
lines, which were obtained from all at present available size measurements, selecting those above 
A > 1.3 cm 2 and performing a fit to (i) only the major axis, (ii) the geometric mean of major 



2 these wavelengths should be least affected by a finite intrinsic source size 



and minor axis, and (iii) the most recent (> 1990) data. The resulting mean slopes in the range 
of (3 = 2.06 ± 0.01 deviate slightly from the slope of = 2 other authors recently have used 
(Bower et al. 2004, Shen et al. 2005 and references therein). Depending on the number of data 
points included in the fit, their weighting scheme and the frequency range over which the fit is 
performed, one obtains slightly different estimates of the intrinsic size of Sgr A*. Since the size 
determination at the shortest wavelengths (A < 3 mm) is still restricted to a circular Gaussian 
component, we determine the de-convolved intrinsic size of Sgr A* using the geometric mean 
size from major and minor axis at all wavelengths. Formally we obtain the following sizes 3 : 
25 ± 7Rs at 7mm, 15 ± 3Rs at 3.4mm, < 65Rs at 2 mm, and 11 =L 6R$ at 1.4mm. 



5. Future Outlook 

Global mm-VLBI performed at the shortest wavelengths (< 3 mm) will allow to study Sgr A* 
(and other AGN like, e.g. M 87) with a resolution of only a few Schwarzschild radii, facilitating 
the direct imaging of the immediate vicinity of super-massive black holes. Global mm-VLBI 
observations of Sgr A* are clearly a great challenge and demand coordinated efforts. In Figure 6 
we show simulations of visibility amplitudes and closure phases for a realistic near future global 
mm-VLBI observation of Sgr A* at 230 GHz. New telescopes like APEX (to be later replaced 
by ALMA) located in the southern hemisphere, and phased interferometers operating as single 
VLBI antennas (PdB, CARMA, SMA) are crucial to provide the necessary sensitivity on the 
longest baselines. When these telescopes become available for mm-VLBI, we may hope to image 
the environment of the SMBH in Sgr A* and detect possible deviations from circular symmetry, 
which are theoretically predicted for a rotating accretion disk black hole system. 
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